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NATIONAL ADVISORY COMMITTEE POE AERONAUTICS 

TECHNICAL MEMORANDUM NO. 1212 

ON THE THEORY OF THE LAVAL NOZZLE* 

By S. V. Falkovich 


In the present paper, the motion of a gas in a plane -parallel 
Laval nozzle in the neighborhood, of the transition from subsonic to 
supersonic velocities is studied. This problem was first consid- 
ered by Meyer (reference 1) who sought to obtain the velocity poten- 
tial in the form of a power series in the coordinates x,y of the 
flow plane. The case of the nozzle with plane surface of transi- 
tion from subsonic to supersonic velocities was further considered 
in a paper by S. A. Christianovich and his coworkers (reference 2). 
For computing the supersonic part adjoining the transition line, 
Christianovich expanded the angle of inclination of the velocity 
and a specific function of the modulus of the velocity in the power 
series, using the velocity potential and the stream function as the 
unknown variables. In a recently published paper, F. I. Frankl 
(reference 3), applying the hodograph method of Chaplygin, under- 
took a detailed investigation of the character of the flow near the 
line of transition from subsonic to supersonic velocities. From 
the results of Tricami's investigation on the theory of differ- 
ential equations of the mixed elliptic-hyperbolic type, Frankl 
introduced as one of the Independent variables in place of the 
modulus of the velocity, a certain specially chosen function of 
this modulus. He thereby succeeded in explaining the character of 
the flow at the point of intersection of the transition line and 
the axis of symmetry (center of the nozzle) and in studying the 
behavior of the stream function in the neighborhood of this point 
by separating out the principal term having, together with its 
derivatives, the maximum value as compared with the corresponding 
corrections. This principal term is represented in Frankl* s paper 
In the form of a linear combination of two hypergeometric func- 
tions. In order to find this linear combination, it is necessary 
to solve a number of boundary problems, which results in a complex 
analysis. 

In the investigation of the flow with which this paper is 
concerned, a second method is applied. This method is based on 
the transformation of the equations of motion to a form that may 
be called canonical for the system of differential equations of 
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the mixed elliptic-hyperbolic type to which the system of equations 
of , the motion of an ideal compressible fluid refers. By studying 
the behavior of the integrals of this system in the neighborhood of 
the parabolic line, the principal term of the solution is easily 
separated out in the form of a polynomial of the third degree. As 
a result, the computation of the transitional part of the nozzle is 
considerably simplified. 


1. Fundamental equations. - The equations of the two-dimensional, 
steady, nonvortical motion of an ideal gas in the absence of friction 
and heat conductivity have the form 


bv _ du 

5x " 5y 


0 


+ S(vp) „ 

ox ay 


0 


w£ JL £ X P 0 
2 + X-l P X-l p 0 


( 1 . 1 ) 

( 1 . 2 ) 


where u and v are the components of the velocity along the x 
and y axes, p is the density, p is the pressure, W = <\/u2 + v2 
is the magnitude of the velocity, x = c p/ c y , Po Pq are 

the density and pressure of the gas at rest. 


Equations (l.l) represent the condition of the absence of vor- 
tices and the equation of continuity. Equation (1.2) is Bernoulli’s 
equation for adiabatic motion for which 



(1.3) 


For the velocity of sound a 



(1.4) 


From equations (1.2), (1.3), and (1.4), the following equation 
is derived: 


P 




(1.5) 
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♦ 


(a Q 2 = Xp 0 /p 0 is the velocity of 
vhich 

dw\ P y 

From equation (l.l), it follows that there exist tvo functions: 
the velocity potential cp(x,y) and the stream function \|/(x,y), 
which are determined by the equations 

dcp = udx + vdy 4}/ = -£- (- vdx + udy) (1.7) 

p 0 


sound in the gas at rest) from 


£0 W. 
p a 2 


( 1 . 6 ) 


In place of the velocity components u and v, the polar 
coordinates, setting u = W cos 9 and v ■ W sin 9, where 9 is 
the angle between the velocity vector and the x axis, are sub- 
stituted. Equations (1.7) are solved for dx and dy, thus 
obtaining 


dx 


cos 9 
W 


p 0 sin 9 


d\|/ 


d y = S3a^a q3+ ^£St®d4, (1.8) 

If x and y as well as W and 9 are considered as func- 
tions of the variables cp and ty, then dx and dy must be 
total differentials, so that the following equations must hold: 

d ( cos e \ _ _ b ( p 0 sin 9 \ 5 / sin 9 \ _ 5 f^O cos 9 \ 

V y \^p W y <5\j/ y W J ckp \jd W J 


In carrying out the differentiation, in taking account of the 
fact that according to equation (1.5) in which p Q /p depends 

only on the magnitude of the velocity W, and in making use of 
equation (1.6), the following equations are obtained: 


sin 9 


5© cos 9 dW 
W Sp 



cos 9 


b9 _ £0 sin 9 
5cp p W 


wf\ dW 

a 2 / 


cos 9 


b9 

S\jir 


sin 9 oW 
W 5\jJ 


5° sin 8 * - £2 11 * 2 . (i - £) 

p Ocp p W \ a z \J ocp 
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By solving these equations for the derivatives 30/3 cp and 
3a/3\J/, 


30 o _ 0 

<5q-> PqW §>}/ 


30 PO a 2 -!? 2 3W ^ 0 
+ p y3 Sq?” 


(1.9) 


This system of differential equations will be of the elliptical 
type if the magnitude of the velocity W is less than the velocity 
of sound and will be of the hyperbolical type for supersonic velocities. 

The new function q is considered instead of velocity W and 
is related to W in the following manner (reference 3): 


"-(i d-io) 

' rr 


Equations (1.9) then assume the form 

- i(n) - 0 t 1 - 11 ) 

where 

. (l -“> 

as a result of (1.10), is a function of the variable t). 

Equations (l.ll) are the fundamental equations for the inves- 
tigation of two-dimensional, nonvortical motion of a gas when the 
velocity of the flow passes from subsonic to supersonic velocity. 

In some cases, It is more convenient in these equations to sub- 
stitute 0 and q as the independent variables and take cp »nd 
\t/ as the required functions. After this transformation, equa- 
tions (l.ll) assume the form 

If + b to) If = ° If ‘ W H = 0 

2. Investigation of variable q. - The variable q deter- 
mined by equation (1.10 ) Is considered in more detail. For 
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computing the integral entering this equation, the square of the 
velocity of sound is 


a = 


k+1 2 


- Sli w 2 

2 W 


( 2 . 1 ) 


i 


In substitutin g the preceding equation in (1.10) 

\2/3 


*1 =12 


1-X 1 

l-X 2 /h 2 ^ / 


h 2 

a* * 


X+l ) 

x-V 


The integration results in 



By expanding equation (2.3) in a series 


*1 


/h(h 2 -l)Y /3 1-X 2 

V 2 ) h 2 -x 2 



( 2 . 2 ) 


(2.3) 


(2.4) 


From equation (2.3), it follows, that T] > 0 for \< 1 and 
tj < 0 for X > 1, that is, in the plane of the variables 9 and 
T), the region lying in the upper half -plane will correspond to the 
region of subsonic velocities and the region lying in the lower 
half -plane will correspond to the supersonic velocities. The line 
of transition from subsonic to supersonic velocity will correspond 
to the line tj = 0, that is, the axis of abscissas. From equa- 
tion (1.10), the value of the velocity W = 0 in the plane 9 ,r\ 
corresponds to an infinitely distant point. For \ > 1, equa- 
tion (2.3) assumes the form 


T) 


(§) 


2/3 


arc 



A 2 -X 2 


- arc tg 



(2.5) 


The characteristics in the plane of the hodograph of the 
velocity for two-dimensional, nonvortical motion of the gas are 
known as epicycloids (fig. 1), the equations of which are (refer- 
ence 4) . 
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9 = C ± h arc tg 



h 2 -* 2 


- arc tg h 



Because for a point transformation characteristics go over 
into characteristics, the following equations of the character- 
istics in the plane of the variables 6 and are found by using 
equation (2.5): 



from which it follows that the characteristics assume the form of 
semicubical parabolas with the cusps on the axis of abscissas 
(fig. 2). 

5. Differential equations of motion of a gas in neighborhood 
of transition line. - The flow in a Laval nozzle near the line of 
transition from the subsonic to the supersonic velocities is con- 
sidered. This line is hereinafter designated the sound line. 

If a straight line perpendicular to the axis of symmetry of 
the nozzle is directed away from the axis, it will intersect the 
streamlines with constantly Increasing curvatures and will there- 
fore encounter particles of the gas having constantly Increasing 
velocity. The sound line will therefore be a curve that is con- 
vex toward the supersonic velocities^- with vertex on the axis of 
symmetry (fig. 3). The point of intersection of the sound line 
with the axis of symmetry is, according to Franhl, denoted as the 
center of the nozzle. 

In the plane of the variables cp and \J/, the region of flow 
is transformed into a strip the width of which is determined by 
the amount of gas flowing through the nozzle (fig. 4) . 

The point of origin of coordinates in the cp,\|/ plane cor- 
responds to the center of the nozzle in the flow plane. 

The determination of the flow reduces to finding two functions 
(cp,\l/) and 0=0(cp,\}/) that satisfy equations (l.ll). Because 
the flow is to be symmetrical with respect to the streamline 
» 0, it is necessary that the required functions satisfy the 
conditions 


Hfhen the streamlines have points of zero curvature , the sound 
line will be a straight line perpendicular to the axis of symmetry; 
this case was considered by S. A. Christianovich (reference 2). 
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Ti(cp,\lO = ti(cp-\j/) s(cp,^) = - 0(cp-\iO n(o,o) = o 

(3.1) 

which are based on equations (2.2). 

The solution of equation (l.ll), in the form of a power series 
in the variables cp and \|/, takes into account equations (3.1) 

2 2 3 2 

q = a^cp + a 2 cp + a 3 \J/ a^cp + agcp\(/ + . . . 

G = bjCf*!' + ^cp 2 ^ + b 3 ^ 3 + b 4 cp^J /+ . . . (3.2) 

from which it follows that if the flow in the neighborhood of the 
origin of coor dina tes is considered, that is, if cp and \|/ are 
assumed to be small magnitudes, the following equations may be 
obtained from equation (3.2): 


rj « o(cp) e = o(cp\j/) = o(\|/) 

l| * °w Ij s o(i) H ■ o(cp) (3 - 3) 

With the use of equation (2.1) and the notations introduced in 
equations (2.2), equation (1.12) for the function b (tj ) may be 
reduced to the form 


b(n) 


fO / h(l-A 2 ) 

P \j(h 2 -X z )r\ 


In accordance with equation (2.4), the following equation is derived: 

b(0) = ( X .l) l/3 = (~f) " (^D l/3 


In taking account of the order of smallness of all terms 
entering equations (l.ll), it is concluded that near the origin 
of coordinates the system of equations (l.ll) may be replaced by 
the following equations: 
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By setting b(0)\j/ , the final result is 

tMt-° <*•*> 

where, for simplicity, the bar over \J/ has been dropped. 

4. Investigation of flow in neighborhood of center of nozzle. - 
It is evident that the functions 


3 2 

e = A 2 cp\{/ - ^ 3 n«Acp-4r^ 2 ( 4 . 1 ) 

where A is an arbitrary constant, are integrals of the system of 
equations (3.4), and satisfy conditions (3.1). 

The significance of the constant A will be explained. From 
the second equation (4.1), q = Acp along the axis of symmetry of 
the nozzle (\) / » 0). Differentiation results in 


A 


_ dq „ dq dW dx 
dqp dW dx dcp 


Furthermore, by using successively equations (1.10) and (2.1) 


*1 „ _ 
dW 



if* / 

* Vitfc 2 -* 2 ) 


Moreover, along the line \J/ =» 0 

dx 1 dW _ 5u 

dcp ~ W dx " 3x 

Bence, for A the following relation is obtained: 


f _h_ [ 1-X 2 du\ 

\ w 2 Vi0> 2 -X 2 ) S 4=o 

7=0 


(X + l) l/5 Au\ 

a* 2 ^ x 'x=0 

7=0 


A = 


(4.2) 
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where to obtain the last result, equation (2.4) was used. 

The value of A is thus proportional to the value of the 
derivative of the velocity at the center of the nozzle. 

It is assumed that 3u/3x >0 so that A will be a negative 
quantity. 

Along the sound line, T] = 0. Hence, according to the second 
equation (4.1), the following equation of the sound line is 
derived: 

(*- 3 ) 

that is, in the plane the sound line will be a parabola. 

From equation (3.4), the differential equation of the char- 
acteristics has the form 

By substituting the value of tj from equation (4.1) 

In the integration of this equation, set 

Acp = | j^y 2 A\| / » x 

The equation then assumes the form 

or 1- g - * 7 

After separating the variables and integrating, the followin g 
equations of the characteristics are obtained: 

x(y+l) 2/3 (2y-l) l/3 = C x(y-l) 2/3 (2y + l) l/3 = C 
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In order to obtain the characteristics passing through the 
origin of coordinates, set C = 0. Thus the variables cp and \J/ 
become 


cp - 


A\J> 2 

2 



(4.4) 


Hence, the characteristics passing through the origin of coor- 
dinates in the cp,\j/ plane are parabolas tangent to each other at 
this point and tangent to the sound line (fig. 4). The origin of 
coordinates will therefore be a singular point of the integrals of 
equations (3.4) determining the flow in the nozzle. 

In considering the character of this singularity, it Is evi- 
dent from figure 4 that the characteristics and the sound line 
divide the neighborhood of the center of the nozzle into six 
regions. It shall be investigated how the neighborhood of the 
center of the nozzle is transformed in the plane of the variables 0 
and T) by the integrals of equation (4.1). By eliminating from 
equation (4.1) the variable cp, the following cubical parabola Is 
used in determining the stream function: 

AV + 3Ati\J/ - 30 » 0 (4.5) 

This equation has one real root if Its discriminant 

6 =» 90 /4 + rj >0 and three real roots If 6 < 0. Because the 

point (cp = 0, \J/ = 0) corresponds in equation (4.1) to the 

point (0 = 0, rj = 0), the equations of the characteristics cor- 

responding to equations (4.1) are in accordance with equation (2.6). 

f 0 


Thus regions I, II, and III of the plane are transformed into the 
same region of the plane 0 ,tj lying between the characteristics 


0» ± §<-T|> 5 / 2 . 


Furthermore, the streamlines \J/ » ± q correspond, as seen 
from equation (4.5), to the straight lines in the 0,ti plane. 
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* 



The transformation of the neighborhood of the nozzle in the 
6,1} plane will thus have the form of the folded surface shown in 
figure 5. The corresponding regions in figures 4 and 5 are denoted 
by the same numbers. 


In order to compute the streamlines in the flow plane, equa- 
tions (1.8) are used in which d\J/ is set equal to 0, after which 
they ass ume the form 


dx = 


cos 6 
W 


d cp 


dy = 


sin 0 
W 


dcp 


By substituting for 0 its value from equations (4.1), the 
magnitude of the velocity W is, according to equation (1.10), a 
function of the variable i} . Thus along the streamlines \J/ = ± q 


dx = rrrrr cos 


wT^T 


. A 2 qcp)dcp dy = ± sin (^~ 


A 2 qcp) dcp 


Integration results in 
Cp 

x = P 7mrr cos 


^ wW 


(^3- - A 2 qcp)dcp 


«P 


AV i a2. 


7 = * \J wItT 8in V3 + A q< ? dcp+ H 


(4.6) 


where H is the width of the nozzle at the critical section. 
In equations (4.6), set according to equations (4.1) 


« A 2 2 
T} = Acp - ^ q 


The computation of the integrals in equation (4.6) reduces, 
evidently, to the computation of the two integrals of the type 



2 

cos A qcp 

w(n) 


dcp 



sin A 2 qcp 

"W^ dcp 
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with the aid of which x and y are expressed as follows: 


A 3 q 3 
I, cos S-SL 


x = Ij_ cos 


3 3 .3 3 

-f- + I 2 .in ±f- 


(< 


y = ± (I 2 sin 


a££ 


i 1 oos St a!) 
(4.7) 


5. Nozzle with surface of weak discontinuity. - The case in 
which weak discontinuities are formed along the Mach lines issuing 
from the center of the nozzle is here considered. For this inves- 
tigation, it is necessary and sufficient that the derivative 
(du/dx) y=0 Possess a discontinuity at the center of the nozzle 

(reference 3) . It is assumed that "both values jjdu/5x) y= qj and 
[(•Vdx) y=0 ] x= are positive. 


x =+0 


From equation (4.1), it is evident that the magnitude A will 
have the value A = A x in the regions VI, V, and IV (fig. 4) and 

the value A = Ag in the region III where, according to equa- 
tion (4.2), A^ < 0 and A 2 < 0. 


From equations (4.1), it is concluded that in the regions VI, 
V, and IV 


2 A l 3 3 

e = Aqcp^ - — \|r 

and for region III 

6 = A 2 2 cp «1/ i-vl' 3 


A, 2 

**1 O 

T| = AjCp - -g-'l' 


Ao 2 

*1 = A 2 cp - — 


(5.1) 


(5.2) 


According to equations (4.4), the equations of the character- 
istics separating the regions IV and V from regions I and II and 
the equations of the characteristics separating regions I and II 
from region III have the forms 


cp = 


A^ 2 

~T~ 



(5.3) 


Substituting the first of these equations in equations (5.1) 
and the second in equations (5.2) 
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A X V 



rp = 


— 


(5.4) 


0 = i- | A 2 V 11 = - A 2 i? 


A oil 2 

CD = - (5.5) 


In order that the flow in the nozzle has no discontinuities, 
it is necessary to determine 6 and t] in regions I and II from 
equations (3.4) in such a manner that the characteristics condi- 
tions, equations (5.4) and (5.5), are satisfied. In order to 
integrate the system, equations (3.4), set 



where f and g are functions to he determined. 

For this substitution, equations (3.4) are transformed into 
a system of ordinary differential equations with the independent 
variable t = <p/ii 2 

2f-2tf ' -g' =0 ff'+3g-2tg' =0 (5.7) 


By the elimination of G' 

g = | [4tf - (f + 4t 2 ) f] (5.8) 


By differentiating equation (5.8) and substituting the result in 
the first equation (5.7), a differential equation of the second 
order for determining f is obtained 


(4t 2 + f) f" + f 


,2 


2tf' + 2f = 0 


(5.9) 


From equations (5.6), (5.4), and (5.5), it follows that the 
boundary conditions for the function f will be 


f 


4 


for t = 


_ A 1 


f = 


Ag for 


t = - -f (5.10) 


In order to integrate equation (5.9), 


it is written in the form 
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- 0 


(The solutions 2tf* - f = 0, that is, f = c*/!T which do not 
satisfy equations (5.10) are eliminated. ) In carrying out the 
quadrature 


f ' + 2t 
2tf ' -f 



f * 


1 f _ 2c l t 
2(t-cj_) t-c^ 


that is, the integration of the linear equation results in 


f » 4cjt - 8c^ 2 


c 2 k /^1 


(5.11) 


t 


The boundary conditions, equations (5.10), which the obtained 
solution equation (5.11) must satisfy, have the form: f = f ^ for 

t = t^ and f = f 2 ^ or t = tg where it is easily seen from 
equation (5.10) that the points (tj, f^) and (t 2 , fg) lie on the 
parabola f = -4t 2 and that t-j_ < 0 and tg > 0. Hence, in order 

to satisfy the boundary conditions, it is necessary from the family 
of parabolas equation (5.11) to choose the parabola passing through 
(ti, f^) and (tg, fg). Upon satisfying these conditions 


C 1 = - 


t l 2 '*~ t l t 2 +t 2 2 
3(t]+ tg) 



lS^-tg) 2 (t 1+ 2tg) 2 (2t x + tg) 2 
27(t 1 +tg)^ 


ft 




It Is necessary that along a streamline the velocity in the 
flow direction should increase monotonically, that is, that q 
should decrease monotonically. Because tj = f\Jr according to 
equation (5.6), f * < 0 must be in the range t^ < t <tg. In 

order to obtain this result, it is necessary that eg < 0. This 
condition Is possible only for 2tj_ + t 2 < 0 and t^ + 2tg > 0 

when in accordance with equations (5.10), the following condition 
is obtained 

A 1 

a !< Ag < — 

for which a flow without discontinuity is possible. 


Translated by S. ReiBS 
National Advisory Committee 
for Aeronautics. 
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